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Abstract 


The amplitudes for virtual Compton scattering off protons are calculated 
within the framework of the diquark model in which protons are viewed as 
being built up by quarks and diquarks. The latter objects are treated as 
quasi-elementary constituents of the proton. Virtual Compton scattering, 
electroproduction of photons and the Bethe-Heitler contamination are dis¬ 
cussed for various kinematical situations. We particularly emphasize the role 
of the electron asymmetry for measuring the relative phases between the vir¬ 
tual Compton and the Bethe-Heitler amplitudes. It is also shown that the 
model is able to describe very well the experimental data for real Compton 
scattering off protons. 
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I. INTRODUCTION 


The understanding of the structure of the nucleon is a fundamental task of particle 
physics. Instead of inclusive reactions where the various nucleon constituents contribute 
incoherently, one may use exclusive reactions to study the structure of the nucleon. The 
simplest exclusive quantities are the electromagnetic form factors of the nucleon, which 
can be measured in elastic electron-proton and/or electron-deuteron scattering. The cross 
section for unpolarized elastic electron-proton scattering is well known up to a momentum 
transfer Q 2 of 33 GeV 2 [!![]. For large momentum transfer (> 8.8 GeV 2 ) only the magnetic 
form factor of the proton has been determined because the contribution of the electric one 
is suppressed in the cross section by 1 /Q 2 . In a recent SLAC experiment || the electric 
and magnetic form factors of the proton have been measured separately in the Q 2 range 
from 1.75 GeV 2 to 8.83 GeV 2 . The data on the neutron form factors is very poor at large 
momentum transfer. 

The next simplest process to probe the structure of the nucleon is Compton scattering 
off protons at large momentum transfers. The cross section for real Compton scattering 
has been measured up to initial photon energies of 6 GeV 0. The process ep —> ep'j 
offers the possibility to study virtual Compton (VC) scattering and its interference with 
the Bethe-Heitler (BH) process. At present there is no experimental data available. In the 
near future, however, CEBAF will offer the possibility to measure the process ep —> epr) 
for incident electron energies smaller than 6 GeV [|3J. The investigation of the high energy 
region requires a new type of electron accelerator ||. We note that the BH contribution to 
the procces ep —> eprf is well known ||. 

The general framework for calculating large momentum transfer processes has been de¬ 
veloped by Brodsky and Lepage 0. Using the QCD factorization theorem for exclusive 
reactions [7].Q, one can write baryon form factors and Compton scattering off baryons as 
multidimensional integrals over a product of distribution amplitudes (DA) and elementary 
amplitudes which represent the scattering of constituents in collinear approximation. A 
DA is the valence Fock-state wave function integrated over transverse momenta k*/ 1 up to 
a scale of order Q 2 . It specifies the distribution of the longitudinal momentum fractions 
the constituents carry inside their parent hadrons. Contributions from higher Fock states 
are suppressed by powers of a s /Q 2 in exclusive processes where a s is the strong coupling 
constant. 

Some information about DAs containing the non-perturbative physics, is obtained from 
QCD sum rules |9|-|TTf and from lattice gauge theory |p~2| . p~3|| . The absolute normaliza¬ 
tion of the nucleon wave function has been estimated: QCD sum rules provide a value of 
In — (5.0 ± 0.3) • 10 _3 GeV 2 for the wave function at the origin of the configuration space 
whereas /at = (2.9 ± 0.6) • lCU 3 GeV 2 has been obtained from lattice gauge theory. Also the 
moments of the nucleon’s DA have been calculated up to order three. These moments are 
used to constrain the first five coefficients in an expansion of the DA over the eigenfunctions 
of the evolution equation which are linear combinations of the Appell polynomials. It is 
obvious that the DAs do not unambiguously follow from the few error-burdened moments 
provided by QCD sum rules or lattice gauge theories (a careful analysis of this problem has 
been performed by Bergmann and Stefanis JL|). Therefore, the DAs to be found in the 
literature, e. g., [pU|,pTj[l4| are to be considered as models. Anyway, the striking feature of 
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these DAs is their strong asymmetry which gives preference to positive helicity u quarks in 
positive helicity protons. The asymmetry of these DAs is the source of theoretical incon¬ 
sistencies in applications of the Brodsky-Lepage model, because it strongly enhances the 
contributions from the end-point regions where a s is large and, hence, the use of perturba¬ 
tion theory is unjustified. Recently it has been shown [ pAj -jlTl that perturbative calculations 
become theoretically self-consistent for momentum transfers as low as a few GeV 2 when the 
transverse separations of the constituents as well as a Sudakov factor, comprising gluonic 
radiative corrections, are taken into account. If one also includes the intrinsic transverse 
momentum dependence of the hadronic wave function perturbative calculations may become 
self-consistent even for momentum transfers as low as 1 to 2 GeV 2 in the case of the pion 
l 8 j and at about 7-9 GeV 2 for the nucleon [[L7]|. On the other hand, the inclusion of the 


transverse degrees of freedom and of the Sudakov factor leads to substantial suppressions of 
the perturbative contributions in the few GeV region which are particularly strong for the 
asymmetric DAs. 

To leading order perturbation theory many processes have been calculated within the 
standard approach of Brodsky and Lepage |7]], among them is the pion’s form factor, the 
magnetic form factor of the nucleon, two-photon annihilation into meson pairs and into 
proton-antiproton as well as Compton scattering off protons. To next to leading order in a s 
only the pion’s form factor and two photon annihilation into meson pairs have been studied, 
ft was shown |R| that the corrections to the lowest order predictions for 77 —> MM become 
sufficiently small (< 25%) only for very large center of mass energies a/s < 10 GeV. Other 
important reactions like elastic proton-proton scattering seem to be beyond feasibility with 
present-day techniques of computing very large numbers of Feynman diagrams. There are 
also unsolved theoretical complications with pinch singularities ||. 

There are two characteristic features of the Brodsky-Lepage model, commonly termed 
the hard scattering approach (HSA): the dimensional counting rules and the conservation 
of hadronic hclicities. The latter feature implies that any helicity flip amplitude is zero and, 
hence, any single spin asymmetry too. The helicity sum rule is a consequence of utilizing 
the collinear approximation and of dealing with (almost) massless quarks which conserve 
their hclicities when interacting with gluons. Whereas the dimensional counting rules are 
in reasonable agreement with experiment, seems the helicity sum rule to be violated by 
20 — 25%. For instance, the polarization parameter P in elastic proton-proton scattering 
is by no means zero. On the contrary, it is rather large (20 — 30%) at moderately large 
values of momentum transfer (px > 2GeV/c ) [{ZCJ]. Even worse, the data show the tendency 
of further rise of the polarization with increasing px- Such peculiar polarization phenomena 
have also been observed in many (moderately) large px inclusive reactions like hyperon 
production in nucleon-nucleon collisions or pp —> [[LLj. The prevailing opinion is that 


also these phenomena cannot be explained in terms of perturbative QCD (see, for example, 
Ref. p2| , |23| ), rather they are generated by an interplay of perturbative and non-perturbative 
physics. 

A model which takes into considerations non-perturbative effects, has been proposed 
by us in a series of papers | 2 ^-^]. In that model baryons are viewed as made of quarks 
and diquarks, the latter being treated as quasi-elementary constituents which partly survive 
medium hard collisions. The composite nature of the diquarks is taken into account by 
diquark form factors. Diquarks are an effective description of correlations in the wave 
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functions and constitute a particular model for non-perturbative effects. The diquark model 
has been applied to a variety of processes and successfully confronted to data. Among these 
applications is a recent study of the nucleon’s electromagnetic form factors [|27]]. In fact, 
predictions are achieved for both the magnetic and the electric form factors. For the latter 
quantity no result is obtained in the pure quark hard scattering model because it requires 
hclicity flips of the nucleon and, in so far, the electric form factor or better the Pauli form 
factor, also represents a polarization effect. In the diquark model helicity flips are generated 
through spin 1 (V) diquarks. The diquark model is designed such that it turns into the 
theoretically well established pure quark picture asymptotically. In so far the pure quark 
picture of Brodsky-Lepage and the diquark model do not oppose each other, they are not 
alternatives but rather complements. 

The purpose of this paper is to discuss virtual Compton scattering off protons and its 
interference with the Bethe-Heitler process where the photon is emitted by the electrons 
(see Fig. 1). The Compton amplitudes are calculated within the framework of the diquark 
model with the parameters taken from previous studies [p6p7|. Real and virtual Compton 


scattering are the two simplest experimentally accessible processes in which the integrals 
over the longitudinally momentum fractions yield imaginary parts. The reason is that there 
are kinematical regions where internal quarks and gluons can go on mass shell. It was 
shown in [281 that the lowest order predictions of perturbative QCD for photon induced 
processes need no resummation of Sudakov corrections. The essential point is that the 
singularities do not pinch as is the case in baryon-baryon scattering. The appearance of 
imaginary parts to leading order in a s is a non-trivial prediction of perturbative QCD. The 
process ep —» epr) with polarized incoming electrons offers the possibility to determine the 
phases of the Compton amplitudes as a function of the center of mass scattering angle, by 
investigating the interference of the virtual Compton with the Bethe-Heitler process. In 
principle, the measurements of other spin observables through the use of polarized targets 
and polarized beams would offer further possibilities of crucially testing the applicability of 
the HSA in the few GeV region. The process ep —> epy has been studied in the framework 
of the pure quark HSA by Farrar and Zhang [|29|. Real Compton scattering has, on the 
other side been investigated by two groups [^9[[30|. Unfortunately, the results obtained by 
the two groups deviate from each other substantially. It is suspected that the discrepancies 
are caused by an improper treatment of propagator singularities in [^9|. The kinematical 
flexibility of ep —> epr) allows to test the predictions with respect to the virtual photon mass, 
to the center of mass energy and to the momentum transfer. 

In Sect. JH] a detailed discussion of the kinematics and the various observables of the 
reaction ep —> epr) is given. A short description of the diquark model for exclusive reactions is 
presented in Sect. Em followed by a brief discussion of the results for real Compton scattering 
(Sect. 0). The predictions for the cross section of VC scattering and for the unpolarized 
cross section of the photon electroproduction ep —> ep'y are discussed in Sect. |Vl This section 
also includes a discussion of the interference between the VC amplitudes and the BH ones. 
The case of polarized electron beams and its physics implications is investigated in Sect. 6. 
The paper terminates with a few concluding remarks (Sect. 7). 
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II. KINEMATICS OF ELECTROPRODUCTION OF PHOTONS 


The process ep —> ep 7 receives contributions from the VC and from the BH processes. 
The kinematics and the helicities are specified in Fig. [L|. We work in the center of mass 
(CM) frame of the final state photon and proton. Neglecting the electron mass the different 
particle momenta can be chosen to be 


fc M = (fc 0 , hi sin a cos fc 0 sin a sin k 0 cosa) 

k ,fJl = (fc 0 — g 0 , k 0 sin a cos /, fc 0 sin a sin /, fc 0 cos a — |p|) 

Qf = (go, 0, 0, |p|) 

q'^ = q' 0 (1, sin 6, 0, cos6>) 

P M = (Po, 0, 0, —|p|) 

P ,fl = (Po, -<?o sin ^ 0, -go cos 9) (2.1) 


in that frame. </> is the angle between the hadronic and leptonic scattering plane, and 6 is 
the scattering angle of the outgoing real photon. The polarization vector e/ of the outgoing 
photon is given by 

e/(//) = ( 0 , — //cos#, — i, //sin 9). ( 2 . 2 ) 

v 2 


go, Po, qh Po an d IpI are related to the Mandelstam variable s = ( p+q ) 2 and to the virtuality 
(g 2 = — Q 2 ) of the exchanged photon by 


?o = 

Po = 


s — Q' 2 — m 2 

VS 

s + Q 2 + m 2 

VS 


P 


% 



s — m 2 
= 2 ^ 
s + m 2 

= 2 ^S 
-g 2 , m 2 ) 


(2.3) 


where m is the nucleon mass. The Mandelstam function A is defined by 


A (x,y,z) = \Jx 2 + y 2 + z 2 — 2 x y — 2 x z — 2 y z. 


(2.4) 


Momentum conservation leads to the following relations 

sin a = ^-J(2k 0 - g 0 ) 2 - |p| 2 , fc 0 cos a = -j—r (g 0 (2 fc 0 - g 0 ) + |p| 2 ) . (2.5) 

21 p | v 21 p | v > 

The other frequently used reference frame is the laboratory frame defined by p = 0. In that 
frame one has the following useful relations: 


Q 2 = 4 k {)L h 0 L sin' 


( 2 . 6 ) 
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where koL and k' 0L are the energies of the inital and final state electrons in the LAB and 
is the electron scattering angle. ^Frorn the theoretical point of view the most significant 
variables are the invariants s, t, Q 2 and the polarization of the virtual photon in the VC 
process 


ST 

0 

1 

to 

0 

to 

1 

p 

2 

(go — 2 ko ) 2 + 

p 

2 


1 + 2 


( k 0L ~ K L ) 2 + Q 2 

Q 2 


tan 



(2.7) 


The unpolarized differential cross section for the reaction ep —> epy can be written as 


d A cT 


ds dQ 2 d(j) dt 


{ep —> ep'f) = 


128 (27t) 4 k% L m 2 A(s, —Q 2 , m?) 4 “ 


1 \ ' \rpBH , rpVC 

7 + Ayy, 


v A 


'x'x 


( 2 . 8 ) 


where T BH and T x c denote the helicity amplitudes for the BH and for the VC contributions 
to the electroproduction of photons, respectively. In the one-photon approximation the latter 
read (see Fig. [I|) 


= A“v»(t , yh"»((:,i') tew x'\%°\px). 


(2,9) 


In the CM frame the two-body s-channel helicity amplitudes := M^ a>a(M,Q 2 ) for VC 
scattering are dehned by 






( 2 . 10 ) 


where ef is the polarization vector of the virtual photon 

ef ( 0 ) = ^ (|p|, 0 , 0 , g 0 ) e“(±l) = +-J= ( 0 , 1 , i, 0 ). ( 2 . 11 ) 

Current conservation leads to a relation between the 0-and 3-components of the hadronic 
current matrix element. The projection of the electron current on the polarization vectors 
is straightforward and leads to the following relation between the VC contribution to the 
photon electroproduction amplitudes and the VC helicity amplitudes 


rjlVC 

\',v\ 


Ana 


Q V 1 ~~ £ ( x/2 - 


\/l + £ + 2 zz \/l — £ e l< ^Mn'x',+ix ~ xf2~£M^ypx 

1 


v/ 2 [ 


\/l + £ — 2 v x/l — e e^M^y-ix) \ 


( 2 . 12 ) 


Due to parity invariance there are only 12 independent s-channel helicity amplitudes con¬ 
tributing to the process 7 *p —> 7 p- These are conveniently denoted by i — 1, 12: 



$5 

- Al + l— i 1 — I 

•e* 

CD 

- M + 1 + i j0+ i 


$6 

= M_i + i x _i 
^ 2 ’ 2 

$10 

= 


$7 

= M 1 1 1 1 

X 2 ,_L 2 

$11 

= M -i+±,o+± 


<f 8 

= M + i + i,i_i 

$12 

= M + i-|,o+i 


(2.13) 
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$1 — $8 correspond to transverse polarization of the virtual photon and $9 — $12 to lon¬ 
gitudinally polarized photons. In the real photon limit the latter four amplitudes vanish. 
In that limit time reversal invariance reduces the number of independent amplitudes even 
further: ($7 = $ 3 , d > 8 = —$ 4 ). The VC contribution to the ep —> epy cross-section can be 
decomposed as follows 


d A a 


ds dQ 2 d(j) dt 
( d&T 


as 


wi¬ 


ve 


x 


dt 


d&L 

+ £ —-b £ COS 2 (p 

dt 


4 ( 27 t ) 2 k^ L m? Q 2 (1 — e) 
d&TT 


dt 


+ \J2e(l + e) cos <f> 


da 


LT 


dt 


(2.14) 


where the partial cross-sections read: 

i) The cross-section for transverse photons (which, at Q 2 = 0, reduces to the unpolarized 
cross section da/dt for real Compton scattering) 


dar c JC 

ir-g' 4 - 

ii) The cross-section for longitudinal photons 


da ,l 
dt 


12 


= cEN 2 


i =9 


iii) The transverse-transverse interference term 
daxT 




iv) The longitudinal-transverse interference term 


(2.15) 


(2.16) 


(2.17) 


‘Gf = ~ Se [$; (4, - * 7 ) + (4 2 + j> 8 ) + 4>;, (4> 3 - <I>,j + $; 2 (4, + 'I',,)]. (2.18) 

The phase space factor c is given by 


c = 


167t(s — m 2 )A(s, — Q 2 , m 2 ) 
The Bethe-Heitler amplitudes read (see Fig. |l|) 

47ra 


(2.19) 


rjlBH 


t L B a H { p y\f\p\) 


( 2 . 20 ) 


where the 7 * p —> p current matrix element is expressed in terms of the magnetic form factor 
of the proton Gm and the Pauli form factor F 2 by: 


(p'\'\j a \p\) = \fknoiu(p , A') ( 7 a G M (t) - y^(p' +p) a F 2 (t)) u(p, A) 


( 2 . 21 ) 
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The electric form factor Ge is related to F- 2 by Ge = Gm — k p (1 + r)F 2 where r = — t/Am 2 . 
The leptonic current L^ H is given by 




+u(k', v') 


ijl To 

Sk'q' 


la jV / 1 

tkq' J 


u(k, v) 


( 2 . 22 ) 


where s a fe = (a + 5 ) 2 and t a b = (a — 6 ) 2 . In the soft photon approximation q' —> 0, the 
first term in ( |2.22| ) is the usual Bremsstrahlung contribution. Putting all together the BH 
contribution to the process ep —» epy may be written in a form similar to that of the 
Rosenbluth cross section for elastic ep scattering: 


d A a 


a 


4 7 r /cql ?« 2 A(s, —Q 2 , m 2 ) t 2 


ds dQ 2 dcj) dt 

x {A(s fcp , s, Q 2 ,4 ?, t) + rG?AJ ( t ^ + ^(gfcp, 5, Q 2 ,0, t) G'lf(t) 


1 + r 

The functions A and i? are given by 


(2.23) 

(2.24) 


A('Sfcp; S-tQ ) ^) ^ 


Sk'q' tkq' 


t 


+ 


fcq' 


Sk'q' 


+ 


-(sfcp - s - Q 2 - 2 


m 


(2.25) 


t 

^kq 


- (s kp + m 2 ) H--— (2 (s fcp - m 2 )(s kp - s - Q 2 ) - u(s + Q 2 ) + m 2 (t + ?n 2 )) , 

' Sk'q' t k q' 


Bi^Skpi s,Q , 1) 


2 4 .\ t ( Sk'q' t tkq' 


2 \ t 


+ 


t 2 Q- 


kq ' 


^k'q' J Sk'q' t k q f 


(2.26) 


For large t the contribution from the electric form factor is suppressed by the factor r as 
compared to that from the magnetic one. The interference term between the BH and the 
VC contribution is a very lengthy expression and we refrain from showing it. 

Similar expressions as (|2.14| ) and (|2.23| ) hold also for polarized cross sections. However, most 
likely only the polarization of the electron beams can be used at high energies. Therefore, 
we refain from discussing polarized cross sections in this article with the only exceptions of 
the electron asymmetry Al obtained by reversing the helicity of the beam electrons, and 
the corresponding proton asymmetry. The discussion of these asymmetries is postponed to 
Sect. 6 . 


III. THE DIQUARK MODEL 

In the hard-scattering model of Brodsky-Lepage 0 the process 7 * p —> 7 p is expressed by 
a convolution of DAs with hard-scattering amplitudes calculated in collinear approximation 
within perturbative QCD. In a collinear situation in which intrinsic transverse momenta are 
neglected and all constituents of a hadron have momenta parallel to each other and parallel 
to the momentum of the parent hadron, one may write the valence Fock state of the proton 
in a covariant fashion (omitting colour indices for convenience) 





















Ip. A) = /sPsOi) Bsti(p, A) + fvfv(x 1 ) Bv(j a + p a /m)j-,u(p, A)/V3 


(3.1) 


w is the spinor of a proton with momentum p and helicity A. The two terms in ( |3.1| ) represent 
conhgurations consisting of a quark and either a spin-isospin zero (S) or a spin-isospin one 
(V) diquark, respectively. The couplings of the diquarks with the quarks in a proton lead 
to the flavour functions 


B s = uS [ U)d ] B v = [uV- V2dV{ UjU y]/V3. (3.2) 


The use of covariant spin wave functions has many technical advantages. For instance the 
calculation of a large set of so-called elementary amplitudes is avoided, one immediately 
projects onto hadronic states. Another advantage is that only hadronic quantities (spinors, 
polarization vectors, and so on) appear. The DA (ps(v)( x i)i where X\ is the momentum frac¬ 
tion carried by the quark, represents the light-cone wave function integrated over transverse 
momentum and is defined in such way that 



dx i v?s, ( vq(ah) 


1 . 


(3.3) 


The constant fs{v) acts as the value of the configuration space wave function at the origin. 
Representative Feynman graphs contributing to the hard-scattering amplitudes for the pro¬ 
cess of interest in this article, are displayed in Fig. [||. The blobs appearing at the gD , 7 gD 
and 77 D vertices symbolize three-, four- and five-point functions. These n-point functions 
are evaluated for point-like diquarks and the results are multiplied with phenomenological 
vertex functions (diquark form factors) which take into account the composite nature of the 
diquarks. Admittedly, that recipe is a rather crude approximation for n > 4. Since the 
contributions from the n-point functions for n > 4 only represent small corrections to the 
final results that recipe is perhaps sufficiently accurate. The perturbative part of the model, 
i.e. the coupling of gluons (and photons) to diquarks follows standard prescriptions, see e. g., 
3T| . The diquark-gluon vertices read (refer to p5|,|26| for notations) 


SgS : ig s tij (pi +_p 2 ) M 

Vgv : -ig s t a ij {g a p{pi+P2)^- g^a [(! + «) p x - kp 2 ] p -g^ [(l + «)p 2 - kpi] q } (3.4) 


where g s = \J\.tux s is the QCD coupling constant, k is the anomalous magnetic moment of 
the vector diquark and t a = X a /2 the Gell-Mann colour matrix. For the coupling of photons 
to diquarks one has to replace g s t a by — a/ 47 raen where a is the fine structure constant and 
eu is the electrical charge of the diquark in units of the elementary charge. The couplings 
DgD are supplemented by appropriate contact terms required by gauge invariance 


qSgS : -2 i e 0 esgstij g, w 

yVgV . i Go6yg s t^j [2 g^ y g a g g^g av g^agfiu) (3-5) 

As we already mentioned the composite nature of the diquarks is taken into account by 
phenomenological vertex functions. Advice for the parameterization of the 3-point func¬ 
tions, ordinary diquark form factors, is obtained from the requirement that asymptotically 
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the diquark model evolves into the hard scattering model of Brodsky-Lepage. This require¬ 
ment fixes the asymptotic behaviour of the form factors. Interpolating smoothly between 
that behaviour and the conventional value of 1 at Q 2 = 0, the form factors are actually 
parametrized as 


Ff(Q 2 ) 


31 

Q 2 s + Q 2 ’ 


4 3 ) (q 2 ) 


( Qv V 
Wv + Q 2 ) ' 


(3.6) 


The asymptotic behaviour of the diquark form factors and the connection to the 
tering model is discussed in more detail in Ref. |EH . 


In accordance with the required asymptotic behaviour the n-point functions for 
parametrized as 


hard scat- 
71 > 4 are 


Fj?\Q 2 ) = asif’W 2 ), 


F ( v\Q 2 ) 


( Qv 

[ av QfTQl 


n —3 

4 3) « 2 )- 


(3.7) 


The constants asy are strength parameters. Indeed, since the diquarks in intermediate 
states are rather far off-shell one has to consider the possibility of diquark excitation and 
break-up. Both these possibilities would likely lead to inelastic reactions. Therefore, we have 
not to consider these possibilities explicitly in our approach but excitation and break-up lead 
to a certain amount of absorption which is taken into account by the strength parameters. 
Actually, for our numerical studies we use |26| , [27 


ips(x 1) = N s XiX 2 exp -b 2 {m 2 q /xi +m 2 s /x 2 ) 
ip v ( Xl ) = N v xxxlil + 5.8 xi - 12.5 Vb r 


1,2 7 ™ 2 




O / \ 


(3.8) 


and the set of parameters 

fs = 73.85 MeV, Q| = 3.22GeV 2 , a s = 0.15, 

f v = 127.7 MeV, Q\ = 1.50 GeV 2 , a v = 0.05, re = 1.39; 


a s = 127t/25 log(Q 2 / Aqcd) is evaluated with A q C d = 200 MeV and restricted to be smaller 
than 0.5. The parameters Qs and Qv, controlling the size of the diquarks, are in agreement 
with the higher-twist effects observed in the structure functions of deep inelastic lepton- 
hadron scattering |52| if these effects are modeled as lepton-diquark clastic scattering. The 
DAs are a kind of harmonic oscillator wave function transformed to the light cone. The 
masses in the exponentials are constituent masses since they enter through a rest frame wave 
function. For the quarks we take 330 MeV whereas for the diquarks a value of 580 MeV 
is used. The oscillator parameter b is taken to be 0.498 GeV -1 . The constant (Ns = 
25.97; Ny = 22.29) are fixed by the requirement (|3.3| ). The more complicated form of the 
DA ipy causes a smaller mean value of x\ than obtained for the DA <ps- The exponentials 
in the DAs provide strong suppressions in the end-point regions. 

According to the above rules the hard scattering amplitudes for the process 7*p —> yp 
are calculated from the set of lowest order graphs of which a few representatives are shown 
in Fig. ^]. In fact, one has to evaluate 32 graphs for each of the diquarks, S and V. For 
comparison, in the pure quark model 366 graphs contribute. The calculation of the 64 
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graphs has been carried out using the SUN version of FORM |33|]. We have checked that our 
amplitudes are gauge invariant with respect to the photons and the gluon. The 7* p —> 7 p 
hclicity amplitudes, decomposed in terms of the various n-point contributions, read 


<*, = 


9 


dm f dxi{fs<p s (yi 


4 $( 5 ’ 3 ) _)_ 2 


<Ps(xi) 


+fv<Pv(yi 


2$P ) -2# 4) + ll# 5) ] <py( Xl )} 


(3.10) 


The factor Cp (= 4/3) is the Casimir operator of the fundamental representation of SU(3) C . 
The n-point contributions generated by the scalar diquarks are explicitely given in the Ap¬ 
pendix |A|. The vector cliquark contributions form very lengthy expressions and we refrain 
form quoting them here; they can be obtained from the authors on request. 


IV. RESULTS FOR REAL COMPTON SCATTERING 


The results of the cliquark model for real Compton scattering obtained from the DAs and 
the parameters in ( |3.8| ) and ( |3.9| ) are shown in Fig. |3] for three different initial photon energies. 
Note that in the very forward and in the very backward regions the transverse momentum of 
the outgoing photon is small and, hence, our model which is based on perturbation theory, 
is not applicable. As compared to the results presented by Kroll, Schiirmann and Schweiger 
25j there are small modifications of minor importance due to the use of covariant spin wave 


functions ( |3 . 1|) and to small changes of the parameters. For the vector cliquark contributions 
the covariant spin wave functions leads to correction terms related to hclicity flips of the 
quarks. Such terms have been neglected in previous work EH . 


Although experimental data are available only at energies which are at the limits of 
applicability of a model based on perturbative QCD, the cliquark model is seen to work 
surprisingly well. The results obtained within the pure quark HSA by Kronfeld and Nizic 
are of similar quality. It is interesting to see that our predictions for the Compton 


cross section do not behave as ~ s~ b at fixed angles and finite energies as the pure quark 
HSA predicts. The reason for the deviations from the scaling law is obvious: the various 
contributions to the cross section exhibit different energy dependences due to the diquark 
form factors. The cliquark model also predicts interesting photon asymmetries and spin 

see the discussion in |^5 

m 


correlation parameters 
proton, of the order of 10%, is predicted 


Even a transverse polarization of the 
This comes about as a consequence of the 


perturbative phases of the amplitudes produced by the propagator poles and of non-zero 
helicity flip amplitudes generated by the vector diquarks. In the pure quark perturbative 
approach transverse asymmetries are zero (for massless quarks). 


V. CROSS SECTIONS FOR VC SCATTERING AND FOR 
ELECTROPRODUCTION OF PHOTONS 

In this section we present the predictions from the cliquark model for VC scattering. We 
select two CM energies, namely s = 5 GeV 2 which is the high energy end of the CEBAF 
accelerator (and is at the same time about the lower limit of the applicability of the cliquark 
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model) and a rather large energy of s = 10 GeV 2 which may become accessible with future 
accelerators like, for instance, ELFE. 

In Fig. [| we show the — * yp cross sections for several values of Q 2 / s. The transverse 

cross section ( [2 .1 5| ) for VC scattering is scaled by s 6 according to dimensional counting. 
The most noticeable fact seen in the transverse cross section is the strong decrease with Q 2 , 
starting from Q 2 = 0, the real Compton scattering case. The effect is particular dramatic 
in the backward direction. It can also be observed from Fig. [4] that, as for real Compton 
scattering, the dimensional counting rule only holds approximately within a factor of 2 be¬ 
tween 5 and 10 GeV 2 . This is a characteristic feature of the diquark model above-mentioned: 
the diquark form factors lead to a transition from a behaviour ~ s 4 (cliquark form factors 
equal unity) to the dimensional counting rule behaviour ~ s~ 6 which is the asymptotic re¬ 
sult. The other three cross sections, oy, &tt and (Jlt, are much smaller than the transverse 
cross section in particular around 90°. Only in the very forward and backward directions 
these cross sections become sizeable but again these regions, say | cos^l > 0.7, are outside 
the hard scattering domain and the application of a model relying on perturbative QCD is 
suspect. In Fig. ^jthe Q 2 dependence of the VC cross sections are shown in a hard scattering 
situation (s = 10 GeV 2 , cos 6 = —0.6) which is accessible at CEBAF. It can be seen that 
the transition from real Compton scattering to large Q 2 VC scattering is non-trivial; a rich 
structure is predicted. A very interesting phenomenon is the change of the slope of the 
transverse cross section. 

The full VC scattering contribution to the cross section for the electroproduction of 
photons is obtained by combining the results shown in Fig. £| or |3| according to Eq. fl2.14| ). 

Next we want to investigate the effect of adding coherently the BH amplitudes to the VC 
ones. To that purpose we plot in Fig. the difference between the full ep — » epy cross section 
and the VC scattering contribution to it divided by the full cross section. The question of 
interest is where are the kinematical regions of small BH contaminations allowing to measure 
the VC process? As can be seen from Fig. dominance of the VC contributions requires 
high energies, small values of | cos0| (the actual values depend strongly on Q 2 and the beam 
energy &ol) and an out-of-plane experiment, i. e., a large azimuthal angle (0 > 50°). For the 
actual CEBAF energy of 6 GeV the VC contribution only dominates in the very backward 
region and for very small values of Q 2 . Outside the regions of VC dominance we expect, 
according to the diquark model, strong BH contaminations (see Fig. j(f). We have however 
not observed pronounced interference phenomena between the VC and the BH contributions 
like the Coulomb-hadronic interference pattern seen in elastic two-body reactions sometimes. 


VI. THE ELECTRON ASYMMETRY 


The regions of strong BH contaminations offer an interesting possibility to measure the 
phases of the VC amplitudes relative to those of the BH amplitudes. As we explain in the 
Appendix some of the internal quarks, diquarks and gluons may go on mass shell. While 
these propagator poles are integrable they lead to phases of the VC amplitudes. Real as 
well as virtual Compton scattering are the simplest reactions in which, to leading order in 
as, such phases appear |2^] and are calculable perturbatively. Thus, it seems to us, that 
measuring these phases is a very interesting check of the hard scattering approach (with and 
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without diquarks). 

As a consequence of the phases of the VC amplitudes, the VC contributions to the 
helicity amplitudes of the electroproduction of photons also obtain non-trivial phases beyond 
the phases due to the azimuthal angle dependence (see Eq. (|2.12|) ). In other words, the 
perturbative phases manifest themselves in the fact that the T matrix is not self-adjoint. For 
the BH process on the other hand, T = obviously holds. As we are going to demonstrate 
information on the absorptive part T — can be obtained from the electron asymmetry 


= o~(+) -o~(~) 

u(+) + a(-) 


( 6 . 1 ) 


where cr(±) are the (differential) cross section for electroproduction of photons with spe¬ 
cific helicity of the incoming electron. Since at CEBAF the electron beam is polarized a 
measurement of the asymmetry seems feasible to us. We emphasize that this possibility 
of measuring the absorptive part of the T matrix by the electron asymmetry does not de¬ 
pend on our model but is quite general and follows from parity and time reversal invariance 
merely. As is well-known a one-particle helicity state transforms under the combined parity 
and time reversal operation as 


|k, A)—> 77(A) |k, -A) (6.2) 

where 77(A) is ±1 depending on A, the spin of the particle and on its internal parity. So the 
combined parity and time reversal operation transform a given helicity state into a state 
with the same momentum but with reversed helicity. If the interaction is invariant under 
the parity and time reversal operation the T-matrix elements for, say, a 2 —> 3 body process 
satisfy the relation 

(k' 1 ,A / 1 ;k',A';k',A'|T|k 1 ,A 1 ;k 2 ,A 2 ) 

= n Vi(K, -A);k', -A';k', —A' \T% 1: -Ap k 2 , -A 2 ;>*. (6.3) 

i 

(For a two-body amplitude one may also use reflection invariance with respect to the scat¬ 
tering plane.) Let us now assume that particle 1 is a spin 1/2 one, say an electron. Then, 
ignoring all kinematical variables, the cross section for particle 1 being in a definite helicity 
state is 


^ . |-^~A , 1 ,A2,Ag;±,A2 1 

{Ai,A'| 

The difference of these cross sections may be written as 


(6.4) 


<r(+) - it(-) = SRe J2 + 

W,A'} * 

' V; M -aJ- 


(6.5) 


If there would be no absorptive part of T, i. e., T — T\ then, according to (|6.3|), the 
difference of the two helicity cross sections and hence the electron asymmetry (|6.1| ), is zero. 
Therefore, the asymmtry measures the non-trivial phase as we claimed above. 
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Specifying this result to the VC contribution, we find for the difference of the two cross 
sections (using the notation of Sect. II and dropping the kinematical factors) 


a(+) -a(-) = - 



sin < 


x 3m [$; ($1 - $ 7 ) + ($2 + $s) + $n ($3 - $5) + ($4 + $ 6 )] (6-6) 

Thus, we see that Aj c measures the imaginary part the longitudinal-transverse interference 
term whereas ctlt (see Eq. ( |2.18|) ) measures its real part. In other words, A\ c measures 
the relative phase between the longitudinal and transverse VC hclicity amplitudes. It turns 
out, however, that the diquark model predicts only very small values for A\ c . The ultimate 
reason is that the longitudinal helicity amplitudes are much smaller than the transverse ones 
for VC scattering (see Fig. [|). So the asymetry is essentially due to the BH-VC interference 
where the transverse VC amplitudes do contribute, thereby producing a sizable asymetry. 
This asymetry is mostly due to the transverse VC amplitudes because, if we switch off the 
longitudinal VC amplitudes, the result changes by only a few percent. 

Another possibility to measure the non-trivial perturbative phase is offered by the proton. 
One may measure the proton asymmetry or the polarization of the (outgoing) proton in 
either the electroproduction of photons or (real or virtual) Compton scattering. For the 
latter two-body processes this asymmetry or polarization is related to the imaginary part 
of the products of hclicity flip and non-flip amplitudes as is well-known. The predictions 
from the diquark model for the proton polarization in the case of real Compton scattering 
is discussed by Kroll, Schiirmann and Schweiger 


I11 Figs. |7| we show the electron asymmetry at a beam energy of 10 GeV and s = 5 GeV 2 
for several values of Q 2 and of the azimuthal angle. Al is generally very small but in 
regions of strong BH contamination it is sometimes spectacularly enhanced. Finally in 
Fig. H we present results for the electron asymmetry for the actual kinematical situation at 
CEBAF (s = 5 GeV 2 k 0L = 10 GeV 2 , Q 2 = 1 GeV 2 ). As can be seen from that figure 
the asymmetry is large in our model for small values of | cosd| (~ 0.2) and values of the 
azimuthal angle around 30°. The magnitude of the effect is very sensitive to details of the 
model and, therefore, should not be taken literally. Despite of this one may take our results 
as an example of what may happen. In so far we believe that is urgent and important to 
explore that phenomenon because it will elucidate strikingly the underlying dynamics of the 
electroproduction of photons in a kinematical situation which can be considered as a (fairly) 
hard scattering region. 


VII. CONCLUDING REMARKS 

We have calculated VC scattering off protons within the framework of the diquark model 
which represents a particular version of the Brodsky-Lepage model appropriate for moder¬ 
ately large momentum transfer. The diquark model combines perturbative QCD with non- 
perturbative elements - the diquarks which represent quark-quark correlations in the proton 
wave function modeled as quasi-elementary constituents. Predictions for the VC-scattering 
cross section and for the ep —> epy cross section are presented for kinematical situations 
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accessible at CEBAF and, perhaps, in the future at an high energy accelerator like ELFE. It 
is also shown that the diquark model predicts a cross section for real Compton scattering in 
fair agreement with the data. We have also elaborated on the BH contamination of the elec¬ 
troproduction of photons which, according to the diquark model, become sizeable for small 
azimuthal angles. The BH contribution offers also the interesting possibility of measuring 
the relative phases between the VC and the BH amplitudes. The phases of the VC ampli¬ 
tudes are a non-trivial phenomenon generated by the fact that some of the internal quarks, 
diquarks and gluons may go on mass shell (note that the ingoing and outgoing quarks and 
diquarks forming the protons are always on mass shell). The electron asymmetry is 
particularly sensitive to the relative phases and can in principle be measured at CEBAF. 
The polarization of the outgoing nucleon in (real as well as in virtual) Compton scattering 
is also a measure for the relative phases between flip and non-flip helicity amplitudes. 

A few words about the self-consistency of our calculations within perturbative QCD are 
in order. As can be seen from ( |A9| ) the running coupling constant a s diverges in the end¬ 
point regions Xi,yi —» 0,1. This is not characteristic of the diquark model, the same happens 
in the pure quark HSA. As a consequence, perturbation theory looses its meaning as a weak- 


coupling expansion. This has previously been pointed out by Isgur and Llewellyn-Smith [34 


in their analysis of form factors. One way out - actually the one we employ - is to freeze the 
running coupling constant once it has reached a certain value (0.5 in the case at hand) ||35|| . 
In addition we make use of DAs which strongly suppress the end-point regions. The freezing 
of a s introduces a new external parameter. The modified HSA proposed by Sterman and 
Li [UHL in which the transverse momenta of the constituents as well as Sudakov corrections 


are taken into account, avoids the introduction of such an external parameter. The Sudakov 
corrections select components of the wave functions with small spatial separations of the 
constituents. With increasing momentum transfer the allowed spatial separations are get¬ 
ting smaller and smaller. Asymptotically, only those components remain which are taken 
into account in the standard HSA of Brodsky and Lepage [7). The numerical effect of the Su¬ 
dakov suppressions is similar to that of freezing in a s . This has been demonstrated recently 
in [HHH for electromagnetic form factors (within the pure HSA). It is claimed that in this 
way the self-consistency of the perturbative contributions can be reestablished for Q larger 
than a few GeV without introducing a freezing parameter. We also note that the transverse 
momentum dependence of the hadronic wave function helps in achieving self-consistency at 
rather low values of momentum transfer [H|. In [^7[ an analogous calculation for the S- 


diquark contributions to the proton form factor has been performed and the authors of |27] 
arrive at the same conclusion. A treatment of the end-point regions in the manner proposed 
by Li and Sterman |T5| diminishes the results for the form factors obtained with the diquark 
model by only 10 to 20%. Since Sudakov corrections mainly depend on colour and not on 
spin a similar behaviour for V diquarks is to be expected. That small suppression can be 
compensated for by adjusting the parameters of the diquark model appropriately. Effects 
of similar magnitude are expected for Compton scattering. We consider these suppressions 
not explicitely but understand them as being absorbed in our parameters. For an effective 
model as the diquark model is, this procedure is sufficient. 
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APPENDIX A: HELICITY AMPLITUDES FOR 7 *qS qS 

The process 7 * q S —> 7 g S only contributes to the hadronic helicity-conserving ampli¬ 
tudes. This implies that 6 of the 12 independent amplitudes < 3y- S ’ n ' ) given in Q2.13Q are zero 


n) = n) = 4>f ’ n) = ’ n) = <hSJ n) = = 0 . 

For dynamical reasons the diquark model yields also 

d4‘ S>) = 0 . 

The remaining amplitudes read in our model: 


>f’ 3) 

4>f ’ 4) = 4 A 


4A^ 3) s 2 


/ XlS 


= —-,- 

ql + ietu(s + Q 2 ) \x 2 y 2 y\u - y 2 Q 2 

(5,4) 1 J (s + Q 2 ) (s + X21) 1 


' uQ- 


2\2 


ym 


-(u + Q 2 ) 


u 


ql + ie gia + i £ alb + i£ 


t + 


U 


$ (5,5) =4A 


(S, 5) 


XiyiU 


y 1 , 


(Al) 

(A2) 


(A3) 


*<*» = 4A? 3) Q 


yiu 


+ 


V\U - y 2 Q 2 u I x 2 y 2 (s + Q 2 ) q\ + ie 


d>< 5 ’ 4 ) = 


4 A?’ 4) 


tQ 2 


1 u + y 2 t 


y 2 u - yiQ 2 u(s + Q 2 ) g 2 a + ie 


$^ s,5) = 


# 3) = 


4Ai 5 ’ 5) 


Q 2 (s + Q 2 + y 2 t) 


y 2 u - y^ 2 xiyiu(s + Q 2 ) 


4 A- 


(S,3) 


Vi U 


+ 


s + Q 2 — 


sQ 2 


24^2(3 + Q 2 ) 


yiu - y 2 Q 2 t \ x 2 y 2 (s + Q 2 ) gf + ie 


s + <5 - 


sQ 2 


(5,4) 


= 4At _ J _1_ 

y 2 u-yiQ 2 \ X]_(s + Q 2 ) 

$(S,5) = 4 A ( t S ’ 5) _ y 2 s 

y2U - yiQ 2 xiyiu(s + Q 2 ) 


+ 


xix 2 (s + Q 2 ) 
x 2 s 1 ) 


(w + y 2 t) _ s_ 

s + Q 2 g\ a + ' xi glb + i£ 


(A4) 


(AS) 


’ 3) = 0 


U2 S 


*(5.4) = 2 a (5,4) 

9 0 (gl + ie)(D 2 1+ is) 

^5,5) = _2 A^ ,5) —-— 


Hi + ie 


*( 5 , 3 ) _ 

^12 — 


*( 5 . 4 ) _ 
^12 — 


4 A 


(5,3) 


yiu 


+ 


ViU - y 2 Q 2 1 1 x 2 y 2 {s + Q 2 ) q\ + 


s + Q 2 


5Q 2 


2 A 


(5,4) 


2s 


y 2 u-yiQ 2 1 xi(s + Q 2 ) 


+ 2 1 + 


V\t 


xix 2 (s + Q 2 ) 
.2 „ „ i 


+ 


2:2s 


s + Q 2 J gj a + i£ xi gl b + i £ 


(A6) 


16 






































^(5,5) _ 2 Aq ’ ^ s s + Q 2 + 2 y 2 t 

12 V2U-yiQ 2 t x^^s + Q 2 ) 


(AT) 


where 


(S» 


-T ~ 

V s e( 

( S,n ) 

V -2 tQ 

-o — 

s 


■<*.{£) Fs\Ql)- 


(AS) 


For Q 2 = 0 the are zero and the remaining amplitudes agree with those quoted by 

Kroll, Schiirmann and Schweiger |25|. The various gluon propagators appearing as argu¬ 
ments in a s and as arguments in the various diquark vertex functions are taken as: 


gl = -x 2 V2 1 g\ = -T 2 {x\y 2 + x 2 yi)t g\ = -xpy\t 


Ql = gl 


Ql = gl + Ql 


Ql = gl + Ql 


(A9) 


(74a, (746, qi and D\ are the momenta of the internal gluon, quark and diquarks which lead 
to the poles within the range of integration. The corresponding virtualities are 


qf = x x s - x 2 Q 2 gl a = xiy 2 s + x 2 yiu - x 2 y 2 Q 2 
D\ = x 2 s - xiQ 2 g\ b = x 2 yxs + x x y 2 u - x x y\Q 2 ■ 


(A10) 


The propagator poles are integrable and do not destroy the validity of the hard scattering 
approach [^]. The poles are handled in the usal way by using the is prescription 


x ± ie 


= P 


x 


=F i 7 rS(x) 


(All) 


In the case that there is only one pole in the region of integration, the integrals over x\ can 
be performed in the following way: 

fl , /(^i,77i,s,t,Q 2 ) 

dx i- 0 , . - 4>d{xi) 


10 


p 


l£ 


_ f 1 , lr J(xi,yi, sp.Q 2 ) (t) D {x i) - f{x^\ T/i, s, t , Q 2 ) (j) D {x S p) ) 


jr 


(A 12) 


+f(xi\y i,s. t . Q 2 ) 4>d(x ?) (pv. j 


1 dx i 


— ITT 


0 p“ 


dp 2 


dx\ 


-i' 


where x\ p} is the zero of the equation p 2 = 0 and the different p 2 are given in Eq. (|A10|) . 
The necessary principal value integral over X\ is then carried out analytically. The different 
integrals and the derivatives of the propagators are listed in Table I. In cases when two 
propagators go on shell at the same time one has to make a partial fractioning before 
applying (|A12|) . 
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TABLES 


TABLE I. Zeros of propagators from Eq. m, partial derivatives and principal value integrals. 
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FIGURES 


FIG. 1. The Bethe-Heitler and the virtual Compton contributions to the process ep —► ep 7 . 
The momenta and the helicities of the various particles appearing in that process are indicated. 


FIG. 2. Feynman graphs contributing to the process "f* p —* 7 p. Graphs with the two photons 
interchanged are not shown. 


FIG. 3. The cross section for real Compton scattering off protons scaled by s 6 vs. cos # for three 
different photon energies. The experimental data are taken from [Q]. 


FIG. 4. The cross section for virtual Compton scattering vs. cos# for several values of Q 2 /s at 
a) s = 5 GeV 2 , b) 10 GeV 2 . Upper left: the transverse cross section scaled by s 6 . Upper right: 
the ratio of the longitudinal over the transverse cross sections. Lower left (right): the ratio of the 
longitudinal (transverse) - transverse interference term over the transverse cross section. 


FIG. 5. The cross section for virtual Compton scattering vs. Q 2 . For notations see Fig. [|. 


FIG. 6 . The difference of the full photon electroproduction cross section and the VC contribu¬ 
tion to it over the full cross section vs. cos 6 for several combinations of values of the beam energy 
ko l, Q 2 and the azimuthal angle <f>. a) s = 5 GeV 2 , b) 10 GeV 2 


FIG. 7. The electron asymmetry Al vs. cos# for several values of Q 2 and of the azimuthal angle. 


FIG. 8. The electron asymmetry at CEBAF. Top: Al vs. (j> for several values of cos#. Bottom: 
Al vs. cos# for several values of cj). 
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